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Abstract 



We analyze the potential of the Next Linear e~^e Collider to study 



^ ■ anomalous quartic vector-boson interactions through the processes e^e 



W^W^ Z and ZZZ. In the framework of SU {2) i®U {\)y chiral Lagrangians, 
we examine all effective operators of order that lead to four-gauge-boson 
interactions but do not induce anomalous trilinear vertices. In our analysis, 
we take into account the decay of the vector bosons to fermions and evalu- 
ate the efficiency in their reconstruction. We obtain the bounds that can be 
placed on the anomalous quartic interactions and we study the strategies to 
distinguish the possible couplings. 



*Email: eboli@ift.unesp.br; ^ Email: concha@axp.ift.unesp.br; ^ Email: mizuka@fma.if.usp.br 



1 



I. INTRODUCTION 



The impressive agreement of the Standard Model (SM) predictions for the fermion- 
vector boson couphngs with the experimental results has been a striking confirmation of 
the S'[/l(2) X f/y(l) gauge structure of the model in that sector However, we still lack 
the same accuracy tests for the structure of the bosonic sector. If the gauge and symmetry 
breaking sectors are invariant under the SUl{2) x f/y(l) gauge group, the structure of the 
triple and quartic vector-boson is completely determined. Thus a detailed study of these 
interactions can either confirm the local gauge invariance of the theory or indicate the 
existence of new physics beyond the SM. 

Presently, we have only started to probe directly the triple gauge-boson couplings at 
the Tevatron ||^,|^ and LEP through the production of pairs of vector bosons. Notwith- 
standing, the constraints on these couplings are still very loose. Future hadron |^ and e"'"e~, 
67, and 77 colliders will provide further information on these couplings and improve 
significantly our knowledge of possible anomalous gauge-boson interactions. 

If the SU {2)l^U{1)y symmetry of the model is to be linearly realized, these studies of the 
triple gauge-boson couplings will be able to furnish information on the gauge-boson four- 
point functions provided that dimension 8 and higher anomalous operators are suppressed. 
This is the case when the breaking of the SU{2)l ® U{1)y symmetry takes place via the 
Higgs mechanism with a relatively light elementary Higgs boson. If, on the other hand, 
no fundamental light Higgs particle is present in the theory, one is led to consider the 
most general effective Lagrangian which employs a nonlinear representation of the broken 
SU{2)l ® U(1)y gauge symmetry. In this case the relation between the structure of the 
three- and four-point functions of the gauge bosons does not hold already at order, 
leaving open the question of the structure of the quartic vector-boson interactions. 

At present the only information on quartic gauge-boson interactions is obtained indirectly 
as they modify the gauge-boson two-point functions at one loop [^. The precise electroweak 
measurements both at low energy and at the Z pole, constrains the quartic anomalous 
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couplings to be smaller than 10~^-10~^ depending on the coupling. 

Direct studies of quartic vector-boson interactions cannot be performed at the present 
colliders since the available centre-of-mass energy is not high enough for multiple vector- 
boson production. This crucial test of the gauge structure of the SM will only be possible 
at the CERN Large Hadron Collider (LHC) through the reaction pp VlVlX PHTU| or at 
the next linear collider (NLC) through the processes e'^e~ VVV ||TTl-|l^, e+e^ FFVV 
e-e- FFVV m, 67 ^ VVF -f-f VV [0], and 77 ^ VVV where V = 
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Z, or 7 and F = e or z/g. 

In this work we analyze in detail the processes e^e~ —>■ W~^W~Z and ZZZ in order 
to assess the potential of the NLC, with and without polarized beams, to study anomalous 
quartic couplings of vector bosons. These reactions will be the most important processes 
to study the quartic gauge couplings at the NLC up to energies of the order of 1 TeV, 
where the processes e~^e~ VVFF start to become important |19[. We work in the 
framework of chiral Lagrangians, and we study all operators that lead genuine quartic 
gauge interactions, i.e. these operators do not give rise to triple gauge-boson vertices, and 
consequently are not bounded by the study of the production of gauge-boson pairs. We 
extend the analysis of Ref. |T^ for the custodial SU{2)c conserving operators taking into 
account realistic cuts and detection efficiencies. Moreover, we also study the non- conserving 
SU{2)c interactions. 



II. THEORETICAL FRAMEWORK 

If the electroweak symmetry breaking is due to a heavy (strongly interacting) Higgs 
boson, which can be effectively removed from the physical low-energy spectrum, or to no 
fundamental Higgs scalar at all, one is led to consider the most general effective Lagrangian 
which employs a nonlinear representation of the broken SU{2)l ® U{1)y gauge symmetry 
pO| . The resulting chiral Lagrangian is a non-renormalizable non-linear a model coupled 
in a gauge-invariant way to the Yang-Mills theory. This model independent approach in- 
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corporates by construction the low-energy theorems that predict the general behavior 
of Goldstone boson amplitudes irrespective of the details of the symmetry breaking mech- 
anism. Notwithstanding, unitarity implies that this low-energy effective theory should be 



valid up to some energy scale smaller than Attv ~ 3 Te V [^] , where new physics would come 
into play. 

To specify the effective Lagrangian one must first fix the symmetry breaking pattern. 
We consider that the system presents a global SU{2)l ® SU{2)r symmetry that is broken 
to SU{2)c- With this choice, the building block of the chiral Lagrangian, in the notation 
of Ref. ||2^, is the dimensionless unimodular matrix field which transforms under 

SU{2)l^ SU{2)r as (2,2): 

= exp (i^^) . (1) 

The (p"" fields are the would-be Goldstone fields and t°- {a = 1, 2, 3) are the Pauli matrices. 
The SU{2)l ® U{1)y covariant derivative of S is defined as 

D^T. ^ d,J: + tg^W^E~tg'J:^B^. (2) 

The lowest-order terms in the derivative expansion of the effective Lagrangian are 



£(2) = ^Tr 
4 



2 .,2 



(D^S)^ (Z^^'S) + f3^g'''- (Tr [TV,])' • (3) 
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where we have introduced the auxiliary quantities T = Sr^S''" and V"^ = (-D^S) S"^ which 
are 5'f/(2)L-covariant and f/(l)y-invariant. Notice that T is not invariant under SU{2)c 
custodial due to the presence of r^. 

The first term in Eq. (|^) is responsible for giving mass to the and Z gauge bosons for 
V = {^/2GF)~^■ The second term violates the custodial SU{2)c symmetry and contributes 
to Ap at tree level, being strongly constrained by the low-energy data. This term can 
be understood as the low-energy remnant of a high-energy custodial symmetry breaking 
physics, which has been integrated out above a certain scale A. Moreover, at the one- 
loop level, this term is also required in order to cancel the divergences in Ap, arising from 



diagrams containing a hypercharge boson in the loop. This subtraction renders Ap finite, 
although dependent on the renormalization scale pO| . 

At the next order in the derivative expansion, D = A, several operators can be written 
down We shall restrict ourselves to those containing genuine quartic vector-boson 

interactions, which are 



4^) = a, [Tr {V,V,)f , (4) 

C't^=a,[TTiV,Vn? , (5) 

C'i^ = a, Tr {VM Tr (TV^) Tr {TV) , (6) 

= «7 Tr iV,yn [Tr iTV^f , (7) 

4? = 1^10 [Tr (TV,) Tr (TK)]' • (8) 

In an arbitrary gauge, these Lagrangian densities lead to quartic vertices involving gauge 
bosons and/or Goldstone bosons. In the unitary gauge, these effective operators give 
rise to anomalous ZZZZ (all operators), W~^W~ ZZ (all operators except C^iq), and 
W^W~W^W~ {df'' and Cf'') interactions. Moreover, the interaction Lagrangians C^q\ 
dj \ and C^iQ violate the SU{2)c custodial symmetry due to the presence of T in their defi- 
nitions. Notice that quartic couplings involving photons remain untouched by the genuinely 
quartic anomalous interactions at the order D = A. The Feynman rules for the quartic 



couplings generated by these operators can be found in the last article of Ref. [^ . 

In chiral perturbation theory, the contribution to the processes e~^e~ W~^W~Z 
and ZZZ arises from the tree level insertion of p'^ operators, as well as from one-loop 
corrections due to the interactions, which renormalize the p"^ operators [^. However, 
the loop corrections to the scattering amplitudes are negligible in comparison to the 
contributions for the range of values of the couplings and center-of-mass energies considered 
in this paper. Therefore, numerically, our analysis is consistent even though we neglected 
the loop corrections and kept only the tree-level p'^ contributions. 
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III. LIMITS ON QUARTIC COUPLINGS 



In order to study the quartic couplings of vector bosons we analyzed the processes 

w+w-z , 



e e 



e e 



zzz 



(9) 
(10) 



which may receive contributions from anomalous WW ZZ and ZZZZ interactions. We in- 
cluded in our calculations all SM and anomalous contributions that lead to these final states. 
Therefore, we consistently considered the effect of all interferences between the anomalous 
and SM amplitudes. The scattering amplitudes were generated using Madgraph in the 
framework of Helas with the anomalous couplings arising from the Lagrangians (§-|D 
being implemented as Fortran routines. Moreover, we include in our calculation the W^s and 
Z's decays taking into account the gauge boson widths, spin structures, and correlations of 
the scattering amphtude. 

We required the visible final state fermions to be in the rapidity region |?7| < 3 and 



separated by Ai? = v^Arp~+~A0^ > 0.7. Furthermore, we also folded in the experimental 
resolution factors associated to the electromagnetic and hadronic calorimetry: 



5E 




0.12 


~¥ 


em 




6E 




0.25 


~e' 


had 


y/E 



0.01 eletromagnetic 
0.02 hadronic . 



(11) 
(12) 



The momentum carried out by neutrinos was obtained using energy-momentum conservation 
after smearing the momenta of final state quarks and charged leptons. As illustration, we 
show in Fig. |T] the effect of the smearing on the reconstructed difermion invariant masses for 
W's and Z's. 

Difermion final states (jj, £^z/, and z/z/; with i = e, fi) were identified as being a W 

or a Z provided their invariant masses after the smearing were in the range respectively. 



0.85Mw , -{Mw + Mz) 



■{Mw + Mz) , l.lSMz 



(13) 



In what follows we present our results for two different center-of-mass energies, 500 
GeV and 1 TeV, assuming an integrated luminosity of 100 fb^^ for both energies. We also 
study the impact of a 80% left-handed and 80% right-handed polarized electron beam while 
keeping the positron beam unpolarized. 

A. e+e" ^ W+W-Z 

We identified W^W' Z events through the topologies 6j, 4j + 2£, 4j + 2z/, and 4j + £ + i/£, 
requiring two difermion systems with invariant masses compatible with the W mass — see 
Eq. (|13D — and one difermion system with an invariant mass consistent with it being a Z . We 
show in Table | the fraction of WW Z events that are reconstructed as WW Z and ZZZ for 
center-of-mass energies of 0.5 and 1 TeV. It is interesting to notice that the reconstruction 
probabilities are basically independent of the e~ polarization. Furthermore, the fraction of 
WW Z events reconstructed as ZZZ generates a background for the study of anomalous 
couplings in ZZZ production. 

The cross section for W^W^ Z [ZZZ) is a quadratic function of the anomalous couplings 

(Ttot = ffsm + E <t + E <no' ' (14) 

i ij 

where asm stands for the SM cross section and 0"°^ (cxano^) is the interference (pure anoma- 
lous) contribution. In Table we present our results after the cuts on rj and AR, but 
before W and Z identification. Therefore, these results should be multiplied by the effi- 
ciencies given in Table |l[ Notice that there are only two independent Lorentz invariant 
structures for the WWZZ vertices at order, which implies that the couplings and ay 
(^4 and ae) give rise to identical contributions to a°^^ and a"no^ in W^W~Z production. 
From this table we can witness that the SM contributions are a slowly varying functions 
of the center-of-mass energy, while the anomalous contributions grow rapidly, as one could 
naively expect. Moreover, the SM background can be efficiently reduced using right-handed 
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electrons as this polarization eliminates almost completely the contribution where the W~ 
couples directly to the e~ fermion line. 

In order to quantify the effect of the new couplings, we defined the statistical significance 
S of the anomalous signal 

S = '""^"^l^^"^' VITe , (15) 
which can be easily evaluated using the parametrization (|T^ with the coefficients given in 



Table fl[ In the above expression e stands for the detection efficiencies presented in Table |. 

Table |ITl| contains the values of the quartic anomalous couplings that lead to an increase 
in the total number of events smaller than 3cr, assuming an integrated luminosity of 100 
fb~^ and that only one anomalous coupling is no n- vanishing. These limits were obtained 
combining events reconstructed as W^W~ Z from all the topologies. It is interesting to 
notice that having right-handed polarized electrons improves the bounds in 20-30% with 
respected to the results for unpolarized beams, while the use of left-handed electrons weakens 
the limits. This result is in agreement with Ref. |T^. Moreover, the bounds improve as 
the center-of-mass energy increases since the anomalous contributions grow with energy. 
In general, more than one anomalous coupling might be non- vanishing. In this case the 
correlation among the anomalous couplings can be easily taken into account using the full 
expression of Eq. (|l^) and Table ||. 

In order to discriminate between the different couplings we studied the kinematical dis- 
tributions of the final gauge bosons. Figure ^ displays the W^W~ invariant mass spectrum 
and the px distribution of the Z in the W^W~ Z production with unpolarized beams at 
^/s = 500 GeV. We plotted in this figure the standard model prediction (dotted line) as well 
as the predictions for = 0.61 (dashed line) and = 0.38 (solid line), which are the values 
that lead to a 3cr signal in the total number of events for unpolarized beams. As we can see, 
the W^W~ invariant mass distribution for presents a larger contribution at small values 
of the W^W~ invariant mass, while gives rise to more events with larger invariant masses. 
In principle we can use this distribution not only to distinguish the anomalous couphngs, 
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but also to increase the sensitivity to the signaL However, this can only be accomplished 
with higher integrated luminosity. On the other hand the pt distributions of the Z are very 
similar in the SM and in presence of the anomalous couplings, being the only difference the 
larger number of events in the latter case. 

B. e+e" ^ ZZZ 

The production cross section for ZZZ final states is much smaller than the one for 
W~^W~Z, and consequently just a few fermionic topologies can be used to identify these 
events. We considered only the final states 6j, 4j + 2i, and 4j + 2z/, and required three 
difermion systems with invariant masses compatible with the Z one according to the pre- 
scription given in Eq. ([T3|). We present in Table |^ the efficiency for the reconstruction of 
the ZZZ final state for the above topologies and center-of-mass energies of 0.5 and 1.0 TeV. 
Analogously to the W~^W~Z case, these efficiencies are independent of the polarization of 
the e~. 

Table |V] contains the values of crgm, cThitJ CTano^ for ZZZ production, taking into 
account the t] and AR cuts, but not the reconstruction efficiencies. At order in chiral 
perturbation theory, all the anomalous interactions are proportional to each other since there 
is only one possible Lorentz structure for the vertex which is multiplied by 04+05+2(0^+0-/+ 
oiq). Therefore, we only present the results for 04, being straightforward the generalization 
to the other cases. From this table we can see that most of the reconstructed ZZZ events 
will be observed in the 6j and 4j + 2z/ topologies. Furthermore, the largest anomalous 
contribution comes from a2no\ being the interference with the SM of the same order of 
the SM contribution but with the opposite sign. Analogously, in W^W~ Z production, the 
anomalous contributions grow substantially with the increase of the center-of-mass energy, 
while the SM cross section decreases slightly. 

We present in Table |V1| the Sa allowed for genuinely quartic couplings that can be 
obtained from the non observation of deviations from the SM in ZZZ production. In our 
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analysis of the 6j topology, we included as background the 6j events coming from W~^W~Z 
that are identified as ZZZ. Despite the reduced number of events in the ZZZ channel, the 
bounds on the quartic couplings are at least a factor of 2 better than the ones drawn from the 
W^W~ Z channel due to the smaller size of the background. On the other hand, contrary 
to the W^W~ Z channel beam polarization does not lead to a substantial improvement on 
the attainable limits. 

We display in Fig. ^ the pseudorapidity and transverse momentum distribution of the 
Z's in unpolarized ZZZ production at -^/i = 500 GeV. As we can see, the anomalous quartic 
interactions leads to more centrally produced Z's (smaller l?]^!) which have a slightly harder 
Pt spectrum. However, the number of reconstructed events is not large enough to allow the 
use of cuts to enhance the anomalous contributions. 



IV. DISCUSSION AND CONCLUSIONS 

W^W~ Z and ZZZ are the best channels for direct study of quartic gauge-boson cou- 
plings in e^e~ colliders with center-of-mass energies up to 1 TeV. At higher energies the 



most e+e^ — W~^W~ff becomes important process [|l4l. We showed in this work that 
the NLC will be able to uncover the existence of anomalous quartic couplings of the order 
(9(10^^) for center-of-mass energies up to 1 TeV and an integrated luminosity of 100 fb~^; 



see Tables |T| and 0. Despite these limits being weaker or of the order of the present indirect 



bounds , the above processes will provide a direct test of the quartic interactions among 
the electroweak gauge bosons. We have also shown that the use of a right-handed polar- 
ized electron beam leads to better limits on the anomalous interactions from the W^W^Z 
production due to the substantial reduction of the SM backgrounds. 

It is also important to devise a strategy to disentangle the anomalous couplings in case 
a departure from the SM prediction is observed. In W~^W~Z production, the analysis of 
the W^W~ invariant mass distribution, see Fig. ^ can be used to distinguish between 
the two possible structures for the WWZZ vertex, one associated to 04^6 and the other 
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related to asj. However, we are still left some two possibilities in both cases. At this 
point it is important to use the information from the ZZZ reaction, because the SU{2)c 
violating interactions leads to a much larger excess of events for the same value of the 
anomalous coupling, due to the coupling structure a4 + Q;5 + 2(a6 + a7 + Q;io)- Therefore, the 
combination of the W^W~ distribution and the total number of events in both reactions 
are a powerful tool to separate the effects of the different anomalous couplings provided 
there is enough statistics. Moreover, the comparison between the W^W~ Z event rates for 
different polarizations can also be used to further distinguish between the couplings a^^Q and 
as J, since the latter are less sensitive to the electron polarization. Finally the anomalous 
coupling ttio has the distinguished characteristics of modifying only the ZZZ production. 
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TABLES 



Topology 


Fwwz% 




6j 


64./52. 


0.8/1.2 


4j + 2£ 


66./55. 


0.8/1.2 


4j + 2i/ 


28./8. 


0.2/0.2 


4 j + £ i/^ 


20./5. 


O./O. 



TABLE I. Fraction of WWZ events that are reconstructed as WWZ and ZZZ for several 
topologies and center-of-mass energies of 0.5/1 TeV. 
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Topology 




O-sm (fb) 


-<t (fb) 


«^ (fb) 


<t (fb) 


«^ (fb) 


«^ (fb) 






5 (1) 


5 (I) 


5 (I) 


5 (I) 


5 (I) 


5 (I) 


6j 


unpol 
pol — 
nnl -1- 


7.41 (8.09) 
13.4 (14.59) 
1 61 (\ 60) 


0.12 (0.12) 

0.0 (0.0) 
22 (0 24) 


2.42 (6.30) 
2.81 (7.38) 
2.0 (5.2) 


0.47 (0.21) 

0.57 (0.) 
34 (0 38") 


4.53 (14.89) 
5.26 (17.36) 
3 77 (\2 A) 


1.18 (-4.71) 
1.37 (-5.52) 
99 f-3 93") 


4j + 21 


unpol 
pol — 
nnl -1- 


0.74 (0.80) 
1.32 (1.43) 
159 fO IB.S"! 


0.01 (0.01) 
0.0 (0.0) 
023 fO 023"! 


0.24 (0.63) 
0.28 (0.72) 
20 fO f)2) 


0.01 (0.02) 
0.05 (0.0) 
03 fO 04") 


0.45 (1.48) 
0.52 (1.73) 
375 (\ 2?,) 


0.12(-0.47) 
0.14 (-0.53) 
01 f-0 39") 


Aj + 2u 


unpol 
pol — 
pol + 


3.06 (2.88) 
5.44 (5.14) 
0.65 (0.58) 


0.01 (0.09) 
0.0 (0.05) 
0.06 (0.15) 


0.94 (3.68) 
1.09 (4.30) 
0.78 (3.08) 


0.22 (0.16) 
0.29 (0.12) 
0.15 (0.18) 


2.09 (6.87) 
2.44 (7.95) 
1.74 (5.7) 


0.66 (-2.73) 
0.76 (-3.21) 
0.55(-2.27) 


4 j +£ui 


unpol 
pol — 
pol + 


5.95 (6.38) 
10.63 (11.6) 
1.28 (1.28) 


0.12 (0.11) 
0.0 (0.0) 
0.2 (0.25) 


1.97 (7.31) 
2.29 (8.60) 
1.64 (6.14) 


0.39 (0.48) 
0.49 (0.43) 
0.26 (0.56) 


3.45 (19.74) 
4.03 (23.12) 
2.89 (16.4) 


0.98 (-3.63) 
1.14(-4.27) 
0.81 (-3.04) 



TABLE II. Values for the standard model, pure anomalous and interference cross sections (see 
Eq. |l^ for the W'^W~ Z production and several center-of-mass energies and polarizations. 
The pol — (pol +) lines correspond to 80% left-handed (right-handed) electron beam polarization. 
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^/s GeV 


e polarization (%) 


"4,6 


"5,7 


500 





(-0.56,0.61) 


(-0.48,0.38) 


500 


-80 


(-0.63,0.63) 


(-0.51,0.40) 


500 


80 


(-0.39,0.49) 


(-0.37,0.28) 


1000 





(-0.38,0.40) 


(-0.26,0.25) 


1000 


-80 


(-0.43,0.43) 


(-0.28,0.28) 


1000 


80 


(-0.27,0.32) 


(-0.20,0.17) 



TABLE III. 3(7 allowed values of the quartic anomalous couplings obtained from the reaction 
e+e" ^ W+WZ. 



Topology 


Fzzz% 


6j 
4j + 21 
Aj + 2v 


59./54. 
62./58. 
25./7. 



TABLE IV. Fraction of ZZZ events that are reconstructed as ZZZ for several topologies and 
center-of-mass energies of 0.5/1 TeV. 
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Topology 




O-sm (fb) 


-<nt (fb) 


«^ (fb) 




^/s/TeV 


5 (!) 


5 (1) 


5 (1) 




unpol 


0.163 (0.145) 


0.169 (0.119) 


2.83 (8.35) 


6? 

J 


pol — 


0.236 (0.21) 


0.224 (0.158) 


3.32 (9.77) 




nn] -1- 


090 I'O 081 


11 I'O 082") 


2 38 (6 98") 




unpol 


0.049 (0.043) 


0.044 (0.034) 


0.846 (2.49) 


4j + 21 


pol — 


0.07 (0.063) 


0.057 (0.047) 


0.985 (2.90) 




pol + 


0.027 (0.024) 


0.03 (0.023) 


0.71 (2.08) 




unpol 


0.204 (0.144) 


0.195 (0.118) 


3.67 (14.2) 


Aj + 2u 


pol — 


0.294 (0.207) 


0.27 (0.156) 


4.29 (16.5) 




pol + 


0.113(0.080) 


0.126(0.077) 


3.08 (11.8) 



TABLE V. Values for the standard model, pure anomalous and interference cross sections (see 
Eq. 14) for the ZZZ production and several center-of-mass energies and e~ polarizations. 



v/^ GeV 


e polarization (%) 


"4,5 


0^6,7,10 


500 





(-0.19,0.25) 


(-0.096,0.12) 


500 


80 


(-0.20,0.26) 


(-0.098,0.13) 


500 


-80 


(-0.18,0.22) 


(-0.088,0.11) 


1000 





(-0.14,0.15) 


(-0.068,0.075) 


1000 


80 


(-0.14,0.16) 


(-0.071,0.079) 


1000 


-80 


(-0.12,0.13) 


(-0.058,0.063) 



TABLE VI. 3a allowed values of the quartic anomalous couplings obtained from the reaction 
e+e" ^ ZZZ. 
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FIG. 1. Reconstructed invariant mass distribution for a jet pair from W and Z decays. The 
full line only includes the effect of the finite width while the dashed line contains also the effect of 
the smearing due to the experimental resolution. 
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FIG. 2. W'^W~ invariant mass and prz distributions for unpolarized W^W~ Z production 
at ^/s = 500 GeV. The dotted line stands for the SM result, while the solid (dashed) line represent 
the case = 0.38 (0:4 = 0.61). 
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FIG. 3. Pseudorapidity and transverse momentum Z distributions in unpolarized ZZZ pro- 
duction at = 500 GeV. The dotted line stands for the SM result, while the solid line represent 
the case 0:4 = 0.24. 
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